Integer quantum Hall effect on a six valley hydrogen-passivated silicon (111) surface 
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We report magneto-transport studies of a two-dimensional electron system formed in an inversion 
layer at the interface between a hydrogen-passivated Si(lll) surface and vacuum. Measurements 
in the integer quantum Hall regime demonstrate the expected sixfold valley degeneracy for these 
surfaces is broken, resulting in an unequal occupation of the six valleys and anisotropy in the 
resistance. We hypothesize the misorientation of Si surface breaks the valley states into three 
unequally spaced pairs, but the observation of odd filling factors, is difficult to reconcile with non- 
interacting electron theory. 

PACS numbers: PAC Nos. 73.40.-c, 73.43.Qt, 71.70.Di 



The silicon field effect transistors (FETs) that are at 
the heart of contemporary microelectronics rely on mo- 
bile electrons or holes confined at the interface between 
Si and a higher bandgap barrier material. This barrier 
in metal oxide silicon (MOS) FETs is SiC>2, an amor- 
phous material which introduces inevitable disorder at 
the Si-SiC>2 interface and limits the carrier mobility in 
these devices. A crystalline interface can be created us- 
ing epitaxial SiGe-Si layers in which mobilities can be 
over an order of magnitude higher than the best MOS- 
FET devices 0, but this technique is limited to the [100] 
oriented surfaces Q- Recently a new technique for Si crys- 
talline interfaces has been demonstrated Q in which a Si 
interface is passivated with a monolayer of hydrogen and 
the barrier material is a vacuum. While the inertness 
and high degree of atomic perfection of these surfaces 
has been known for some time 0, 0,13] > the development 
of high mobility electronic devices on H-Si enables the 
exploration of two-dimensional (2D) physics of novel Si 
surface orientations and may one day allow quantum de- 
vices to be engineered at the atomic scale using surface 
manipulation techniques 0. 

We report here the first detailed magneto-transport 
studies of a 2D electron system (2DES) at a H-Si(lll) 
surface gated through a vacuum barrier. Electron mobil- 
ities are an order of mag nitude higher (24, 000cm 2 /Vs) 
than Si(lll) MOSFETs, enabling the observations of 
the integer quantum Hall effect (IQHE). In the effec- 
tive mass approximation the ground state for a 2DES 
on the Si(lll) surface is six-fold degenerate with each Si 
conduction band valley contributing an equal number of 
carriers, (Fig. Id) each with anisotropic in-plane masses, 
m x = 0.19m o and m y — 0.67m o Q (m is the mass of 
the free electron). However, measurements of Shubnikov- 
de Haas (SdH) oscillations in Si(lll) MOSFETs have 
shown conflicting valley degeneracies of two[l, [ij [l(| 
and sixpj], [l2| along with isotropic resistivities for both. 
Subsequent proposals 11 , 13[ have tried to explain these 
anomalous observations, but to date conclusive experi- 
mental results are still lacking. 

The high mobility 2DES can be created by contact 
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FIG. 1: (a) Schematic cross-section of a H-Si(lll) substrate 
contact bonded to a SOI substrate. A p + layer in the SOI 
defines the gate, where blue arrows depict the electric field. A 
2DES is formed at the H-Si(lll) surface within an encapsu- 
lated cavity, (b) The H-Si(lll) substrate has four n + contacts 
numbered accordingly. Tilted magnetic fields are applied in 
the x-z plane, (c) A 1 /im x 1 /im AFM image of atomic 
steps on a H-Si(lll) surface in relation to the crystal direc- 
tions and the contacts of the device, (d) The projection of the 
six valleys for the Si(lll) surface with pairs of valleys labeled 
A, B and C. 

bonding two individual Si substrates Q (Fig. la). One is 
the H-Si(lll) substrate (float zone, p-type, p ~ 10 Sl-cm) 
which has four phosphorous contacts forming a 1-mm- 
wide square with sides oriented parallel to the [110] and 
[112] crystallographic directions (Fig. lb and lc). The 
second is a silicon-on-insulator (SOI) substrate which 
acts as the remote gate, where an electric field can be 
controlled within an etched cavity. The Si(lll) surface 
is H-passivated by immersion in an ammonium fluoride 
solution. The two substrates are then bonded in vacuum 
(^10 -6 Torr), which allows the remote gate to induce 
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TABLE I: Energy gaps for various filling factors, where 8 — 0° 
corresponds to a magnetic field normal to the 2DES. The 
subscripts in A(i2,34) represent the measurement orientation. 

electrons on the H-Si(lll) surface and protects the air 
sensitive surface inside the cavity. 

Figure 2a shows traces of the Hall (Rxy) and the lon- 
gitudinal resistance (R xx ) as a function of perpendicu- 
lar magnetic field, B±, for a constant electron density, 
n = 6.5 x 10 11 cm" 2 at T = 300 mK. We have obtained 
similar data in two devices over a wide range of n. We 
present data for a single device at a density where IQHE 
features are most apparent. There are two orientations 



in which R xx is measured; R 



(12,34) 



and i?(i3.24), where 



the first and second subscripts represent the current and 
voltage contacts respectively (Fig. lb). Anisotropy in 
R xx is observed, but both orientations exhibit minima 
that occur at integer values of filling factor, v — nh/eB±. 
Similar behavior (within 3%) is also observed for R(43,2i) 
and 11(42,31)- These minima are directly correlated with 
the observed plateaus in R xy , a trademark of the QHE. 
Filling factors less than six are observed at v = 5, 4, 
3, and 2, where v = 5 is apparent through d 2 R xx / ' dB\. 
Figure 2b shows characteristic behavior at low B fields: 
Rxy/B± is dependent on B±; both Rn2,34,) an< ^ -^(13,24) 
display a positive magneto-resistance (MR); finally, while 
more apparent in Jim, 34), both orientations display SdH 
oscillations with minima occurring approximately every 
eighth Landau level (LL) filling factor, Av ~ 8. 

Classical MR calculations[14] show that in a multi- 
component electron system in which electrons have dif- 
fering anisotropic in-plane masses, R xy is non-linear with 
respect to B± , and R xx should display a positive MR. In 
Figure 2b we plot R xy /B± and R XX (B) overlaid with a 
classical MR model in which all six valleys are equally 
occupied, with n and the mean scattering time, t, as free 
parameters. The values which provided the best fit to 
R xy are within 5% of observed values: r = 3.5 ps and 
n ~ 6.8 x 10 11 cm~ 2 . This model agrees well with the 



observed behavior in R xy and the positive MR in R xx . 
However, in order to explain the observed anisotropy in 
R xx , an unequal population of electrons among the six 
valleys is required. 

The strong R xx minimum at v = 2, suggests a 
gap between the six valleys. We obtain gap energies, 
A„, from the temperature dependence of R XXl R xx ~ 
exp(— A^/2/cbT). These measurements show that the en- 
ergy gap, A2, is insensitive to small changes in B± and in- 
plane fields, Bn (Table 1). It is by far the largest gap (~ 7 
K) and provides strong evidence that the six valleys are 
indeed broken asymmetrically following a 4-2 valley con- 
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FIG. 2: (a) R xy (green), R, 
versus B± for the device at T 



(12,34) 



(blue) and -R(i3,24) (red) 
300 mK and n = 6.5 x 
Plateaus in R xy are labeled by horizontal bars 
and the numbers indicate their corresponding filling factor v. 
The inset is a schematic energy level diagram depicting the 
three energy scales affecting T>(E) at B± = 1 T. (b) The 
same data from (a) is plotted within < B± < 2 T. Hall 
resistance is plotted as R xy /B± and the divergence is removed 
for B± < 0.1 T. Open symbols (R xy /B±) and closed symbols 
(R xx ) correspond to calculations having all six valleys equally 
occupied. The inset depicts T>(E) for the 7 K model at B — 0. 



figuration. For simplicity we will assume a model which is 
composed of only one type of splitting which is indepen- 
dent of B: two valleys (A) are 7 K lower in energy than 
the remaining four valleys (B and C) in the ground state 
(inset of Fig. 2b). The inset of Fig. 2a illustrates the ef- 
fect B± has on the the density of states, T>(E), for this '7 
K model'. Quantization of the electron's orbital motion 
creates LLs which are separated by the cyclotron gap, 
Ec = heB±/rriT>, where mp = [m x m y Y^ 2 — 0.358m o is 
the density of states mass for Si(lll). Each LL then splits 
into two energy levels due to spin. This energy separa- 
tion is the Zeeman energy, Ez = g* hbB±/ cos6, where 
g* is the effective g- factor and is the Bohr magneton. 

The validity of this 7 K model can be tested with mea- 
surements in R xx versus tilted B fields. Rotating the 
device some angle, 9, away from B± changes Ez with 
respect to Ec, and eventually energy levels of different 
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FIG. 3: Data and calculations represented at n = 6.5 x 
10 11 cm -2 and T = 300 mK. For the color scales repre- 
sented in (a), (c) and, (d), a maximum in d 2 R xx /dB 2 ± or 
dEp/dn corresponds to a minimum in R xx . (a) A surface 
plot of d 2 _R( 12j 34)/d^ 2 versus l/cosS for filling factors within 
5 < v < 11. At 8 C , the R xx minimum at v = 8 disappears 
while a minimum at v — 10 is at a peak, (b) Schematic energy 
level diagram following the 7 K model predicting the two-fold 
energy levels from the A valleys (spin up and spin down) 
crossing one another at 8 C . (c) A plot of d 2 R( 12 ,3A) / dB\ 
versus v and l/cos#. The overlaid schematic fan diagram 
(dashed lines) represents expectations from the 7 K model 
with Az/ ~ 8. (d) A plot of a calculated dEF / dn versus v and 
1 / cos 9 using the 7 K model. The same fan diagram used in 
(c) is overlaid on top of calculations. 

spins will overlap with one another. The angles at which 
these coincidences occur can be observed through the dis- 
appearance of minima in R xx . One such coincidence is 
shown in Fig. 3a, where a minimum in d 2 i?(i2,34) jdB\ 
corresponds to the disappearance of the R xx minimum 
at v = 8 and a concomitant maximum occurs at v = 10. 
Similar behavior is also observed for d 2 i?( 13 _ 2 4)/dBj_. 
The second derivatives of R xx are plotted in order to 
amplify structure and to eliminate the contribution from 
the positive MR at low B fields. Figure 3b depicts a 
schematic energy level diagram following the 7 K model 
and using B± values centered at v — 8. The 7 K model 



identifies this observed behavior with a crossing between 
the two-fold degenerate A valleys in upper spin level of 
the lowest LL (N — 0) and the lowest spin level of the 
second LL (N = 1). This crossing occurs when Ez = Ec- 
Assuming the effective mass equals m-p = 0.358to o , we 
obtain an enhanced g* — 3.3 (bare value g — 2). When 
these parameters are used, the 7 K model can explain 
the observed strengthening of the v = 6 minimum with 
increasing 9 prior to the level crossing and, although not 
shown, explains why v — 12 is never observed. 

The 7 K model with m-p = 0.358to o , implies ~ 50% 
of the electrons are occupying the lowest two valleys at 
small B fields. If only these contributed to the low B 
SdH oscillations, a periodicity of Av ~ 8 would be ex- 
pected if the degeneracy due to spin and valley pairs 
with opposite k remain unresolved. Figure 3c shows 
the remainder of the coincidence measurements, where 
d 2 R(i2.34)/dB^ is plotted versus l/cos# (x-axis) and v 
(y-axis). The plot is a compilation of -R(i2,34) measure- 
ments oriented at 45 different #'s and evidence of up to 
~ 5 level crossings is apparent. A schematic energy fan 
diagram in accordance with expectations from the 7 K 
model, Ais — 8 and g* — 3.3, is overlaid on top of the 
data, and the intersections of the levels coincides with 
minima in d 2 R{\2,34) I ' dB\ or maxima in R xx . Figure 3d 
shows calculations derived from the 7 K model of the 
derivative of the Fermi energy with respect to the den- 
sity, dEp/dn, in the same range plotted in Figure 3c. 
The same fan diagram used in Figure 3c is also overlaid 
on top of the calculations. By applying different Dingle 
temperatures, 0.6 K and 2 K, to the carriers originating 
from the two-fold and four-fold valleys respectively, the 
7 K model is able to reproduce the unresolved spin de- 
generacy and the absence of oscillations from the higher 
energy valleys seen in the data. Differing scattering times 
are commonly observed in two component electron sys- 
tems with different densities [151]. 

A likely origin for breaking the valley degeneracy in 
our 2DES is the miscut of the wafer. This miscut angle, 
tjj, can be derived from the width of the atomic steps on a 
H-Si(lll) surface. Figure lb shows an atomic force mi- 
croscope (AFM) image of a H-Si(lll) surface obtained 
from an adjacent piece of the Si(lll) wafer used to fabri- 
cate the device. From the AFM image we observe steps 
running at an angle, <fi ~ 15 ± 5°, off from the [112] axis 
with an average width of ~ 90 nm. This corresponds 
to a miscut of tp ~ 0.2° and is approximately about the 
[110] axis. In the effective mass approximation, this mis- 
cut increases the normal mass, m z , of the two valleys 
(A) aligned with the [112] direction by ~ 0.5% relative 
to the remaining four valleys (B and C). As a result, the 
six-fold symmetry is broken into a 4-2 valley configura- 
tion similar to our 7 K model, but with an energy gap 
of only A4_2 ~ 2 K. Additionally, the four valleys (B 
and C) occupied at higher energies are split due to the 
off axis rotation of <f). This energy gap is calculated to 
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temperatures are also observed in other high mobility 2D 
systems 3, lj| and have been attributed to correlated 
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FIG. 4: R xx minima at v = 3, 4, and 6 are divided by their 
respective R xx values at 8 = 0° and plotted within —15° < 
8 < 25°. Data represents measurements at T = 300 mK and 
n = 6.5 x 10 11 cm -2 . 

be much smaller, A4 ~ 0.2 K, which is within our exper- 
imental uncertainties, but potentially explains the pres- 
ence of v — 4. It is possible that the discrepancy in the 
magnitude of A2 is due to a many body enhancement. 
However, these predicted energy gaps do not include pos- 
sible effects the atomic steps at the surface may have on 
the confinement of the electronic wavefunction[l6(. 

The simple 7 K model does break down at high B, 
specifically for the observed minima at v = 5, 4, and 3. 
Figure 4 displays the behavior of the R xx minima for 
v = 6,4, and 3 versus small 9. In agreement with the 7 K 
model, the v = 6 minimum strengthens as ~ 1/ cos 9 for 
small positive and negative tilts about 9 = 0°. However, 
the R xx minima for v = 4 are observed to have a linear 
response to small tilts about B±. For surface electrons 
which have non-zero off diagonal terms normal to the 
surface in the effective mass tensor, it is possible to have 
confinement energies depending linearly on 9 17j . The 
magnitude of this effect is oc sin 9 cos 9 and has a peak 
value of ~ 1.5 K for v = 4. This agrees well with the 
observed linear response around small 9 for v — 4. 

The observation of v = 3 and 5 indicates that the 
in-plane symmetry of ± k|| between opposite valleys is 
broken on our H-Si(lll) surface. We also observe v = 3 
to be the only state which exhibits anisotropy in A be- 
tween the two resistance orientations (Table 1). Figure 
4 shows that this anisotropy is also dependent on the 
measurement orientation with respect to £>|| . This split- 
ting cannot be explained by uniform strain or misori- 
entation effects in the effective mass approximation. It 
also cannot be easily interpreted using theories developed 
for valley splittings on [100] Si surfaces [j[ Il6j|. where the 
confinement potential couples valleys with minima in the 
direction normal to the surface. It is possible that disor- 
der on the surface couples the opposite k|| valleys in our 
samples. However, anisotropies of R xx developing at low 



electron phases. The physics may also be similar to the 
lifting of the sublattice symmetry recently observed in 
graphene[2ol|, or to skyrmion states in multivalley elec- 
tron systems [2ll |22| . 

We note that our sample fabrication process is not 
yet optimized and that prospects for substantially im- 
proved mobilities are good. Our results imply that the 
valley splitting on the Si(lll) surfaces can be tuned by 
wafer orientation and in-plane magnetic fields. Thus, H- 
passivated Si surfaces can be prepared with the maximum 
number of valley degrees of freedom in which to explore 
potentially new many-body quantum phenomena, such 
as exotic fractional quantum Hall states [23j|. and 'valley 
singlet' states which are potentially relevant for topolog- 
ical quantum computing [24]. Lastly, the relative ease of 
incorporating surface engineering within the fabrication 
process may enable future experiments in which high mo- 
bility electrons are coupled to molecules or atoms on an 
atomically clean semiconductor surface. 
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